Noise suppression in reproduction of single flux quantum pulses 
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We present the analysis of the mean switching time and its standard deviation of an overdamped 
Josephson junction, driven by a direct current and a single flux quantum (SFQ) pulse. The per- 
formed analysis allows to find the optimal value of the bias current of the clock generator, responsible 
for the shape of SFQ pulse, which minimizes noise-induced switching errors. 



The investigation of the fluctuational characteristics of 
Josephson junctions (JJ) is very important due to their 
applications as logic devices Existence of fluctua- 

tions in Josephson junctions leads, for example, to limit- 
ing lifetimes of the information unit, recorded in Joseph- 
son memory cell, to random switching of logic gates and 
to spread of arrival time of signals in Josephson transmis- 
sion lines (see, for example, and Q). In the present 
paper we consider the dynamics of a short overdamped 
JJ driven by a direct current and SFQ pulse. An SFQ 
pulse is supposed to be born by another JJ named a clock 
generator. We present the analysis of the mean switch- 
ing time (MST) and its standard deviation (SD) of the 
overdamped JJ versus bias current of the clock genera- 
tor, that is responsible for the shape of SFQ pulse. The 
performed investigation allows to find the optimal value 
of the bias current of the clock generator, which mini- 
mizes switching errors induced by noise. In addition we 
test the limits of applicability of the formula for the SD 
recently derived by Semenov and Inamdar 0. 

It is well-known, that in the frame of Resistively- 
Shunted- Junction model Q a point Josephson junction 
in the limit of a small capacitance (high damping), driven 
by current I with fluctuations taken into account is well 
described by the Langevin equation: 



_id<p(t) du(ip) 
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here u(tp) = 1 — cos ip — i(t)ip is the dimensionless 
potential profile, ip is the difference in the phases 
of order parameter on opposite sides of the junc- 
tion, i(t) = io + f(t) = I/I c , I c - critical current of JJ, 
ip (t) = Ip/I c , Ip is the random component of the cur- 
rent, uj c = 2eR]\[I c /h is the characteristic frequency 
of the junction. In the case when only thermal fluc- 
tuations are taken into account, the random current 
may be represented by white Gaussian noise: (ip(t)} = 

2-7 

0, {ip{t)i F {t + t)) = — S(t), where 7 = 2ekT/M c is 

the dimensionless intensity of fluctuations, e is the elec- 
tron charge, k is the Boltzmann constant, T is the tem- 
perature, and h is the Planck constant. 

Initially, a current going across the junction is smaller 
than the critical one, io < 1- A single flux quantum 
pulse arriving from a clock generator (or simply from an- 
other J J) switches the junction to the resistive state, that 



leads to generation of new SFQ pulse. This pulse will be 
generated not immediately, but at the later time which is 
called the switching time. Since due to noise the moment 
of pulse generation is a random quantity, let us investi- 
gate its mean and standard deviation. As a clock pulse 
we choose a voltage pulse, shape of which can be obtained 
by solution of equation Q without any fluctuations 

for i > 1: f(t) = A , "Xl ■ 7 + 1 - a), 

where a is the current going across the junction utilized 
as a clock generator, u — \fa? — 1 is the oscillation fre- 
quency, ip — 7r + arctan(— u>), A is the signal magnitude. 
In the inset of Fig. 1 the form of the current pulse 
= io + f(t) is presented for the case of the equal 
maximal value of current i(t) = 2.9. 
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Fig. 1: The MST versus current a for 7 = and different 
bias currents: io = 0.9, A = 1 - line with squares; io — 0.8, 
A = 1.05 - solid line; io = 0.7, A — 1.1 - dashed line. Inset: 
the form of current pulse i(t) for A = 1.05, io = 0.8 - solid 
line, A = 1, io = 0.9 - dashed line. 



By the definition the first and the second 

moments of switching time take the following form: 



(t) 



tw(t)dt, (t 2 ) 



t 2 w(t)dt, (2) 



where w(t) — dP(t)/dt, P(t) is the probability to find 
tp whithin the interval (— ir, it) and the standard devi- 

In the following 



ation is, as usual, a = y (t 2 ) — (t) 2 . 
t and a will be determined both via direct computer 
simulations of Eq. Q and by numerical solution of the 



corresponding Fokkcr-Planck equation using the Crank- 
Nicolson scheme. In all figures both r and a are normal- 
ized to 1/lUc- Note, that for presently used technological 
processes l/ui c is of the order of 1 ps (see @>0)- 

In Fig. 1 the MST is presented versus current a in the 
case of zero noise intensity for three different values of 
bias current io, but constant total current. MST has a 
minimum as a function of current a, which for larger bias 
current io becomes more broad and deep: minimal MST 
increases of about 1.2 times for bias current io — 0.7 in 
comparison with io = 0.9. 
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Fig. 2: The MST and SD versus current a for different 
noise intensities: 7 = 0.1 - dashed line, 7 = 0.05 - solid 
line, 7 = 0.01 - line with diamonds, io = 0.9, A = 1. 
Inset: io = 0.7, A = 1.1, 7 = 0.001 - line with squares, 
7 = 0.01 - solid line, 7 = 0.1 - dashed line. 



In Fig. 2 MST and SD versus current a are presented 
for different noise intensities for io = 0.9, A = 1. As it 
is seen, these characteristics have minima as functions of 
a, and the minimum of r is shifted to larger values of 
current a, and is more broad and flat than the minimum 
of a. One can see that for noise intensity 7 — 0.05 the 
most optimal range a is located around 1.3-1.4, because 
here the minimum of a is reached, and MST is near the 
minimum and does not almost depend on the noise inten- 
sity (curves t for 7 = 0.1 — 0.01 around a = 1.4 actually 
coincide). Increasing of MST minimum with increasing 
of noise intensity is due to the noise delayed decay ef- 
fect that has been studied in 0, in connection with 
Josephson junctions. More clearly this effect appears, if 
one takes smaller bias current i — 0.7 and changes noise 
intensity in broader interval (7 = 0.1, 0.01 and 0.001). As 
it is seen in the inset of Fig. 2 the location of MST min- 
imum and also its value significantly depends on noise 
intensity: minimal MST for 7 = 0.1 is about 2 times 
larger than the minimal MST for 7 — 0.001. Therefore, 
high-r c devices must be redesigned with account of this 
effect to get maximal performance. 

In Fig. 3 MST and SD are presented versus generator 
current a for different values of bias current io and signal 
magnitude A, but constant total current i + 2A = 2.9. 



It is seen, that for larger bias current io switching occurs 
faster and is less random: minimum of a is wider for 
i = 0.9 and A = 1 than for i = 0.7 and A = 1.1. It 
is necessary to point, that minimal SD value does not 
depend on bias current, and is determined by the noise 
intensity only (see also Fig. 4). 
T, a 
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Fig. 3: The MST and SD versus generator current a (7 = 
0.02): io = 0.9, A = 1 - line with squares; io = 0.8, 
A = 1.05 - solid line; io = 0.7, A — 1.1 - dashed line. 

In Fig. 4 the SD versus generator current a is presented 
for io = 0.9 and noise intensities 7 = 0.1; 0.01; 0.001. One 
can see, that for 7 < 0.01 minimal SD value scales as the 
square root of noise intensity as for the case of a step- wise 
signal |3]-0- The asymptotic expression of SD, derived 
in for the case 7 C 1 (that is somewhat more exact 
than one, derived in Q) has the form: 
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The asymptotic values of SD, given by ©, are presented 
in Fig. 4 as dashed straight lines: for 7 < 0.01 they 
are close to the minimum of <r, the disagreement does 
not exceed 20%. In paper it has been demonstrated, 
that the jitter can be minimized by the use of the sharp 
pulses. As it follows from the performed analysis (see Fig. 
4), even for soft SFQ pulses the jitter can be suppressed 
down to the level of sharp driving case if the parameters 
are properly optimized. 

Fig. 5 shows the results of computer simulations for 
values of bias current and noise intensity that are typ- 
ical for real RSFQ circuits (A = 1, 7 = 0.001, i = 
0.5; 0.7; 0.9). As it is seen, for bias current equals 0.5 and 
0.7 the minimum of standard deviation is not observed 
in the considered range of parameters. This could be a 
reason why this effect had not been observed before. 
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Fig. 4: The SD versus current a for io = 0.9, A — 1 
and different values of noise intensity, from top to bottom 
7 = O.f; 0. Of; 0.001: circles - computer simulations, solid 
curves - formula dashed straight lines - formula 0. 
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Fig. 5: The SD versus current a for different values of bias 
current for 7 = 0.001. Circles, diamonds and triangles - 
computer simulations for io = 0.5; 0.7; 0.9, respectively. 
Dashed, solid and long-dashed curves - formula IgJ. 



Very recently a universal, but approximate formula for 
the jitter a has been derived in |4J. Since this formula 
has neither been compared with the results of computer 
simulations nor investigated versus bias current of the 
clock generator we have performed a certain analysis to 
test its validity for the considered task. In our notations 
the formula for the jitter by Semenov and Inamdar 
may be presented in the form: 
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junction and 7 is the noise intensity. The formula has 
been derived with the assumption that fluctuations adia- 
batically follow the MST in working point that results in 
a jitter of switching time. In spite of the simplicity of 
it gives surprisingly good coincidence with the results of 
computer simulations, see Fig.s 4 and 5, especially for 
relatively large noise intensity 7 = 0.1 and smaller bias 
currents 0.5 and 0.7. For bias current i$ = 0.9 and noise 
intensity 7 = 0.001 the error may be up to 50 — 100% 
and the minimum of a versus a observed in the present 
paper is not reproduced. 

We have considered fluctuational dynamics of a short 
JJ driven by a direct current and SFQ pulse. It has been 
demonstrated that both mean switching time and its 
standard deviation have minima as functions of the bias 
current of a clock generator, responsible for the shape 
of SFQ pulse. Therefore, by proper choice of the bias 
current of the junctions, both the response time and the 
jitter of the rapid single flux quantum logic devices can 
be minimized. 
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where r is the MST, io < 1 is the bias current of the 
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